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Abstract 

It is shown that it is possible to establish sum rules that must be 
satisfied at the nodes and extrema of the eigenstates of confining po- 
tentials which are functions of a single variable. At any boundstate 
energy the Schroedinger equation has two linearly independent solu- 
tions one of which is normalisable while the other is not. In the domain 
after the last node of a boundstate eigenfunction the unnormalisable 
linearly independent solution has a simple form which enables the 
construction of functions analogous to Green's functions that lead to 
certain sum rules. One set of sum rules give conditions that must 
be satisfied at the nodes and extrema of the boundstate eigenfunc- 
tions of confining potentials. Another sum rule establishes a relation 
between an integral involving an eigenfunction in the domain after 
the last node and a sum involving all the eigenvalues and eigenstates. 
Such sum rules may be useful in the study of properties of confining 
potentials. The exactly solvable cases of the particle in a box and 
the simple harmonic oscillator are used to illustrate the procedure. 
The relations between one of the sum rules and two-particle densities 
and a construction based on Super symmetric Quantum Mechanics are 
discussed. 

PACS: 02.30. Gp, 03.65.-w, 11.30.Pb 
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1 Introduction 



There is a well defined procedure for constructing Green's functions for de- 
scribing solutions to linear second order differential equations with inho- 
mogeneous terms (Morse and Feshbach 1953 [1]). This procedure may be 
employed to study solutions to the Schroedinger equation in one dimension. 
For problems with spherical symmetry partial wave decomposition effectively 
reduces the three-dimensional Schroedinger equation to a radial equation in 
r-space and hence the techniques for constucting Green's functions in one- 
dimension are applicable. The Green's functions may be used to establish 
trace formulae in which the integrals over the Green's function may be related 
to sums involving the eigenvalues of the homogeneous differential equation 
(Berry 1986 [2], Sukumar 1990 [3], Voros 2000 [4]). Such trace formulae are 
very useful in checking the accuracy of the numerically computed eigenval- 
ues. The recent interest in the real spectra of non-hermitian Hamiltonians 
exhibiting PT symmetry (Bender and Boettcher 1998 [5], Mezincescu 2000 
[6], Bender and Wong 2001 [7]) has initiated accurate numerical computation 
of eigenvalues of PT symmetric Hamiltonians and the trace formulae have 
proved to be useful. 

In this work we consider the two linearly independent solutions to the 
Schroedinger equation at an eigenenergy and show that it is possible to con- 
stuct functions which are suitable for studying various sums over eigenstates 
in the domain outside the last node of a chosen eigenfunction. In section 2 
it is shown that it is possible to construct new sum rules involving all the 
eigenstates and eigenvalues. In section 3 the examples of a particle in a box 
and the simple harmonic oscillator are used to illustrate the sum rules. The 
relations between the sum rules, two particle densities and Super Symmetric 
Quantum Mechanics are discussed in section 4. Units in which h — 1 and 
the mass ji — \ are used throughout the paper so that ^- — 1. 
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2 Nodes, Extrema and Sumrules 



We develop a formalism in this section which would be suitable for applica- 
tions to either the radial Schroedinger equation for a specific partial wave in 
the domain [0, oo] or the full line [—00, +00]. We consider the solutions to 

= ( V - E >) *J (!) 

satisfying the boundary conditions at the lower and upper end points of the 
domain at x and x 1 

Lt r ^j ->0, Lt r ^ xl Vj -> (2) 

for the normalised eigenstates corresponding to the eigenvalues Ej. The 
state with j = 1 corresponds to the groundstate with no nodes and the state 
tyj has (j — I) nodes. Let the nodes of the eigenstate be at r = R and 
let the outermost node be at r = R . Then \P n has no nodes in the domain 
i?o < r < x\. The second linearly independent solution at E n is then 
given by 

*" (r) = *" (r) IW(ij* (3) 

where R is a constant which may be chosen according to some appropriate 
requirement. If we choose R > Rq then in the domain R < r < x\ there 
can be no infinities arising from the denominator inside the integral and 
is well defined in this domain. The Wronskian relation 

■to - \j> = 1 (A) 

ar ar 

shows that Lt r ^ Ro 7^ but has a finite value determined by the derivative 
of \l/ n at r = R . 

The differential equation satisfied by the eigenstates may be used to rep- 
resent the Wronskian between the states j and k (j 7^ k) in terms of the 
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overlap integrals between the different orthonormal eigenstates in the form 



^ k (y) (y) dy = I ^ k (y) ^ (y) dy = ^— — - -L- (5) 

x J xi \ h k ~ Ej) 

where the dot denotes a derivative with respect to r. If we now define 



where the sum is over a complete set of eigenstates excluding the state n, 
then using the equality in eq. (5) it can be established that 

-*n(r)) G(r,f)= 5>,-(f) f *„ (y) (y) dy . (7) 



The Green's functions G(r,f;E) considered in the usual textbooks (Morse 
and Feshbach[l], for example) are constructed at energies E which are not 
one of the eigenenergies E n . In contrast the function G considered here is 
constructed with E = E n . Using the completeness relation satisfied by the 
eigenstates 

X>,-(f) *i(y) = ^-V) ~ *n(r) * n (y) (8) 

it can be shown that 

*n(r)^ - *n(r)^ G (r, f) = -* n (f) (q (r - f) jf +9 (f- - r) ^ 2 n dy 

(9) 

where 6*(z) is the unit step function which vanishes when z < and has value 
1 when z > 0. 



2.1 Sumrules at nodes of \I/ n 

Various interesting relations follow from the differential equation (9). If we 
choose r = R Q where Rq is a node of the state \l/ n at which ty n (R ) = then 
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we get the relation 

(10) 

In particular setting f = Rq in eq. (10) leads to 

^ {En - Ej) 1 > 

Squaring both sides of eq. (10), integrating over the variable f in the full 
range [xo, xi] and using the orthonormality of the states \Efj it may be shown 
that 

At any node of any eigenstate the rest of the eigenstates must fulfill the 
conditions implied by eqs. (11) and (12). 



2.2 Sumrules at extrema of ^ n 

Another special case of of eq. (9) arises when r — R 1 where R\ is an ex- 
tremum of \l/„ which satisfies (R\) = 0. Eq. (9) simplifies to 



(E n - Ej) 



Ri J xq 



(13) 

which when differentiated with respect to f and evaluated at f = R\ leads 
to the identity 

Using the completeness relation of the eigenstates the above relation may 
also be given in the form 



Squaring both sides of eq. (13), integarting over the variable f in the full 
range [xoj^i] and using the orthonormality of the states \&j it can be shown 
that 

^ ^ 2 (Ri) 1 r Ri r xi 

Swf^-W. •■■*/. (16) 

At any extremum of any eigenstate the derivative of all other eigenstates 
must satisfy the conditions implied by eqs. (15) and (16). 



2.3 Integral realtions valid outside the last node of ^ n 

The differential equation (9) satisfied by the function G defined by eq. (6) is 
a first order differential equation which can be brought to the form 

d G (r, f) * n (f 



dr m n (r) 



( e {r ~ ° r * 2jz r) l ^ idz ) (i7) 



and can be integrated from an arbitrary point r 2 to get the relation 

Gjryr) G(r 2 ,f) 
V n (r) ^ n (r 2 ) 



(f) f ^ ( Hy - f) r -° (f - y) n 

J r2 n V J y J xo / 

(18) 



It is also possible to establish a differential equation for G(r,r). Using 
eqs. (5) and (6) and the limiting value Lt z ^ 9 (z) = 1/2 it can be established 
that 

^ ^ 2 (r) 

S(r) ee G(r,r) = £(^4^ (19) 



satisfies 



d S 



n d~rn = U Kdy ' J X0 Kdy] ■ m 



This equation can be integrated from any point r 2 to give 
S(r) = 



li sM+ *" w (II(f^-£^))- (21) 
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Different choices of f and r 2 in eq. (18) lead to different integral relations. 
If we set f = r 2 in eq. (18) then the resulting expression can be rearranged 
to give 



G(r,r 2 ) G(r 2 ,r 2 ) f r dy 



m n (r) tf n (r 2 ) *2( r2 ) J r2 W n (y) 



9(y-r 2 ) -9(r 2 -y) H 9 2 n (z)dz. 

(22) 



By interchanging the labels r and r 2 another relation like the one given above 
may be derived and by addition of the two relations we can show that 



E 



^(r 2 )V f r> dy 



^E n -E 3 \y n (r) * n (r 2 )J J r< V*(y) (23) 



where r< (r>) is the smaller (larger) of (r, r 2 ). The integrand in eq. (23) is 
free of singularities in the domain of integration if both r and r 2 are greater 
than the last node R of ^ n . 

Various integral relations follow from eq. (23). For example multiplying 
eq. (23) by ^ (r) (r 2 ), integrating over both the variables from R to x 1 
and using the notation 

Ajk = ! 1 *j (y) (y) dy (24) 
and noting that when r = R is a node of \l/ n eq. (5) gives 

^ - HIT^ • ' * " ' (25) 

we can establish that 

a- E - *i (A.) E ^ - - /: f ^ f »; w * 



./\ — it — i .i 



f Xl dr r f Xl 
(26) 
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A special case of the above relation arises if we consider the groundstate with 
n — 1 for which R = x and for all the eigenstates \I>j (^°) = ^- We thus 
get the sum rule 

gCT-f ™*r wt)J?^ dz (27) 

which expresses the inverses of the separation of the eigenvalues of a confining 
potential from the groundstate eigenvalue in terms of an integral over the 
nodeless groundstate eigenfunction. 

The main results derived in this paper are the relations expressed in eqs. 
(9), (11), (12), (15), (16), (23), (26) and (27). In the following sections 
exactly solvable examples will be used to illustrate the sum rules derived in 
this section. 

3 Examples of sumrules at nodes and extrema 
3.1 Particle in a Box 

In this section we consider the example of a free particle confined in a box 
with infinite walls at xo = and x\ = it. The normalised eigenfunctions and 
eigenenergies are given by 

*j(r) = ^1 sin jr , E j = f , j = 1,2... . (28) 

There is a node of the eigenfunction \P n at R = it (n — 1) jn. We first 
examine the sum 

v ' tt ^ (n 2 - ? 2 ) 

which can be simplified using partial fractions, addition formulae for trigono- 
metric functions and standard sums over sine functions (Gradshteyn and 
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Ryzhik 1965 [8]) to the form 

sin nr sin nf ( 1 

G(r, r) = h r cot nr + r cot nr — 7T cot nr> (30) 

7rn \ 2n 

where r> is the larger of (r, f) . Using 



lh *l(<l)d q ^ — (31) 
, 7T n 

and 

^n(i?o) = y|ncos(n- 1)tt = (— l) n_1 n \j~ (32) 
it is simple to show that for r = R eq. (30) becomes 

^'^S^"^^"^"^) (33) 

thereby verifying eq. (10). For the choice f = R = n — n/n eqs. (29) and 
(33) can be used to give 

= <«) 

verifying eq. (11). Eq. (33) can be squared and integrated over f to show 
that 

sin 2 jR n 2 n — 1 



V^ = ^ (35) 

which is the sum rule arising from eq. (12) in this case. 
We next examine 

G(r,r) = - }_^J 2 — ( 36 ) 

i+n J 

which can be evaluated by taking the derivative of the relation in eq. (30) 
with respect to r. There is an extremum of \I/ n at R\ — n — ir/ (2n). Hence 

sin tit 

G (R u f) = (— i?i esc nRi + ir esc nRi 9 (R 1 - f)) . (37) 

7T 



Using 



I 



Ri p 

*ldy = § (38) 



'0 * 

it can be shown that 
-9(f-R 1 )j\ 2 n dy + 6(R 1 -f)J\ 2 n dy= (-^- + 9 (R ± - f)) (39) 

which together with eq. (37) verifies eq. (13) at the last extremum of 
at R 1 . Differentiation of eq. (37) with respect to f, evaluation at the point 
f — Ri and use of the completeness relation leads to 

2 ^-^ . 2 COS 2 ji?l 2 >r-^ .9n ,An\ 

^ j£ (™ 2 - J 2 ) * j 

It is possible to prove this directly by starting from the equality in eq. (30) 
for r = f and show that for Ri = n — ir/ (2n) 

\g{r u r 1 ) = E^-rf- = "4^» 

«fd 2 n( A, J 2 cos 2jgi 3 



which leads to the relation 

+ sin2ji?i ) = - 1 = - §in2ni?i (42) 

thereby verifying eq. (15). By squaring eq. (37) and integrating over f it 
can also be shown that 

^ j 2 cos 2 jRi tt 2 2n- 1 

^(rf-f) 2 16 n 2 ^ 
which is the sum rule arising from eq. (16) in the present case. 

We next examine eq. (21) which in this example becomes 

G (r r) G (r 2 r 2 ) f r dy ( f n f y \ 

A = - - ' - = / — tt- I / sin 2 nz dz — / sin 2 nz dz | . 

Kir) nfa) J r 2 sin 2 y \J y J J 

(44) 
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Using eq. (30) it can be shown that 

a 2r — n 2r 2 — 7T . . 

A = cotnr — cotw 2 (45) 

2n 2n y J 

in agreement with the direct evaluation of the integral on the right hand side 
of eq. (44). 

We next examine eq. (23) which in this example gives the relation 
G(r,r) G(r 2 ,r 2 ) G(r,r 2 ) f r > _dy_ _ _tt_ 

n (r) n (r 2 ) *n (r) * n (r 2 ) ~ 2 J r< sin 2 y ~ 2n [C ° tm> C ° tnr<) 

(46) 

Using eq. (30) to express the various terms on the left hand side of eq. (46) 
it is easy to check that the sum of the terms on the left hand side yields the 
expression on the right hand side of the equation. 

The triple integral on the right hand side of eq. (27) for this example can 
be evaluated to give 

2 r . 2 J r d y r . 2 , 3 

— / sm x ax / — / sm z az — - (47) 
Jo Jx sm y J y 4 

and using partial fractions it may be shown that 

oo 

Ztt^—I (48) 

thus verifying eq. (27). 



3=2 v J 



3.2 Simple Harmonic Oscillator 

We consider an oscillator potential V = x 2 in the range [—00,00] correspond- 
ing to a frequency cu — 2. The oscillator length parameter equals 1 in the 
units we have used in this paper. The energy levels and the eigenf unctions 
are given by 

E j+1 = (2j + 1) , = (^j ' ^/^exp (-x 2 /2) Hj (x) , j = 0, 1, 2, ... 

(49) 
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where Hj(x) are Hermite polynomials which satisfy 

^ = 2j H^x) , H 2j (0) = (-Y M , H 23+ M = . (50) 

We examine the sum rules arising from the choice n = 2 which corresponds 
to the first excited state with has a single node at x — 0. All the antisym- 
metric states with even values of j vanish at x = and the symmetric states 
corresponding to odd values of j have limiting values at x = given by 

*^ (0) = \H {wm) (ir) 2 '^ ' 1 "" (51) 

where the first two factors on the right hand side arise from the normalisation 
integrals of the harmonic oscillator eigenf unctions (Pauling and Wilson 1935 
[9]) and the last factor arises from the values of the even order Hermite 
polynomials at x = (Abramowitz and Stegun 1965 [10]). Hence 



* 2 (Q) 



v 




(52) 



which verifies eq. (11) for n = 2. 



We next examine 

,2 



y> ^ (0) jl (2fc)! 1 

j^{E 2 -E j f ^ V 167T ^(2^!) 2 {2k -if 



Lt z ^ (zarcsinz + (l - z 2 ) 1/2 ^j 



8" <»> 



The normalised eigenfunction for n = 2 given by 

1/4 



/4\ 1 

q> 2 ( x ) = I - J x exp (-x 2 /2) (54) 
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can be used to show that 



J nO poo 

-j— - / Vjdy / *\dz = 

m(o) J-oo Jo 



4 («) 



1(0) 

which when considered together with eq. (53) verifies eq. (12) for the n = 2 
first excited state of the simple harmonic oscillator. 

To examine the sum rule arising from extrema of eigenfunctions we con- 
sider the groundstate n — 1 which has an extremum at x — 0. For all the 
symmetric states corresponding to all odd values of j the derivative of the 
eigenfunction at x = vanishes and for the antisymmetric states correspond- 
ing to even values of j the derivative at x = is given by 

^-wM 1 ?)''^ 1 ' 2 '-- (56> 

Using the values of the eigenfunctions and their derivatives at x = given 
by eqs. (51) and (56) we can show that 




- = -*?(0) (57) 



7T 

thereby verifying eq. (15) for the groundstate of the oscillator. 



We next consider 

y> g (0) 1 (2k)\ 

f^iEi-E,) 2 V4vr ^(2^!) 2 (2* + l) 




l A^T 

— Lt z ^i (arcsin^;) = J— . (58) 
It can be shown that for the normalised groundstate eigenfunction 

*i (*) = (- J exp (-x 2 /2) , ^ y *?dy ^ * 2 dz = i v^F (59) 
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which when taken together with eq. (58) verifies the sum rule given by eq. 
(16) for the n — 1 groundstate of the oscillator. 

4 Discussion 

In this paper sum rules which must be satisfied at the nodes and extrema 
of boundstate eigenfunctions of confining potentials have been established. 
The sum rules in eqs. (11), (12), (15) and (16) have been verified for the 
case of a particle confined in a box and also explicitly for the case of a simple 
harmonic oscillator in the ground or first excited states. However the sum 
rules are valid for all states of the oscillator and for all confining potentials. 
When scattering states are present the expressions have to be modified by the 
addition of an integral to take account of the contribution from the scattering 
states to the sum over the contribution from the discrete states. 

We have shown that in the domain after the last node of an eigenfunction 
the feature that the inverse of the eigenfunction is singularity free may be 
used to establish a variety of relations between the values of all the other 
eigenfunctions in this region and integrals involving the nodeless eigenfunc- 
tion. We have illustrated the sum rules in eqs. (23) and (27) for the case 
of a particle in a box for which the sums and integrals converge and can be 
carried out analytically. For the harmonic oscillator the sum and integral in 
eq. (27) do not converge. 

The antisymmetric wavefunction for two non-interacting identical fermions 
moving in the same single particle potential V such that one of them is in 
the state \l/ n and the other in tyj is given by 

^■(ri,r 2 ) = J\ (tf„(n) V s (r 2 ) - * n (r 2 ) (n)) . (60) 
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The relationship in eq. (23) may also be given in terms of $ nj in the form 



which sheds an interesting light on the sum rule in terms of joint probabilty 
density of two-particle states. 

Supersymmetric Quantum Mechanics may be used to interpret the inte- 
gral on the right hand side of eq. (26). If we consider a potential V which 
is identical to V in the region outside the last node of \I/ n at Rq but has 
an infinite wall at the last node, then the boundstate solutions in V must 
vanish at r = R and as r — > x\. The groundstate energy of V must be 
Ei = E n because \&„ goes to zero at r = R and as r — > x\ but has no nodes 
inbetween. \P n is the groundstate eigenfunction of V but has to be renor- 
malised to 1 in the region [R , xi\. Let the other boundstate eigenvalues of V 
satisfying boundstate boundary conditions at R$ and x\ be Ej, j = 2,3, .. . 
A supersymmetric partner to the potential V constructed by the elimination 
of its groundstate at E\ = E n is 



which is free of singularities for r > Rq. This construction which is based 
on the methods of Supersymmetric Quantum Mechanics (Sukumar 1985 [11]) 
guarantees that the boundstate spectrum of V\ is identical to that of V except 
for missing the groundstate of V at E\ (i.e) V\ has spectrum Ej,j = 2, 3, .... 
It may be shown that a solution at the energy E n in V\ is $ = l/^ n - From 
this solution two other solutions which satisfy boundary conditions at Rq and 
x\ can be constructed in the form 




(61) 



Vi = V - — 2 ln* n (r), r > R , 



(62) 




(63) 
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with the Wronskian 

W = $! ^ $ 2 - $ 2 ^ $! = jf 1 ^(y) dy . (64) 

These solutions may be used to construct a Green's function for the potential 
V given by 

G 1 (r,r>r) = MlMl , Lt r ^ Q G - 0, Lt^G - (65) 

and (5i(r, f) = Gi(f, r). The trace of this Green's function (Sukumar 1990 
[3]) is related to the spectrum of V\ by 

/. Gi(r,r)dr = - / ^ / ^dy / = £ — . (66) 

Using eqs. (60), (61), (64) and (66) it may be shown that 

i r Xi r Xi r Xi 

V- — / / ^ 2 nj (r 1 ,r 2 )dr 1 dr 2 =W / Gi(r,r)dr 

= ( r\ 2 n (y)dy] 1 — 

(67) 

expressing the trace of the Green's function for the Supersymmetric partner 
potential V\ with the energy spectrum (Ej,j = 2, 3, ..) in terms of two-particle 
densities in the potential V with the energy spectrum (Ej,j = 1, 2, ..n, ...). 

We conclude by reiterating that the sum rules expressed in eqs. (11), (12), 
(15) and (16) must be satisfied at all the nodes and extrema of the boundstate 
eigenfunctions of confining potentials in 1-dimension. Also the sum rule in 
eq. (23) for confining potentials is a key result derived in this paper. As 
noted before it is possible to extend the sum rule to potentials which have 
scattering states by the addition of an additional integral to include the 
contribution from the scattering states in addition to the contribution from 
the discrete states which are included in eq. (23). We have focussed attention 
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on confining potentials because of the existence of exactly solvable problems 
for which the sum rules can be explicitly verified. We have verified the sum 
rules for two exactly solvable confining potentials. We have interpreted one 
of the sum rules using the notion of two particle densities and established 
a connection with the trace of the Green's function of a Super Symmetric 
partner in Super Symmetric Quantum Mechanics. 
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